Universal properties of the Higgs resonance in (2+l)-dimensional U(l) critical systems 
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We present spectral functions for the magnitude squared of the order parameter in the scaling limit 
of the two-dimensional superfiuid to Mott insulator quantum phase transition at constant density, 
which has emergent particle-hole symmetry and Lorentz invariance. The universal functions for the 
superfiuid, Mott insulator, and normal liquid phases reveal a Higgs resonance which is relatively 
sharp and is followed by a damped oscillation (in the first two phases only) before saturating to 
the quantum critical plateau. In order to understand the counter-intuitive Higgs resonance in the 
insulating and normal phases, we invoke a picture of a scale-dependent Mexican hat. Our results are 
derived from analytically continued correlation functions obtained from path-integral Monte Carlo 
simulations of the Bose-Hubbard model. 
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Field theories of a complex scalar order parameter, VP, 
can have two types of collective excitations. The first 
one originates from fluctuations of the phase of iff and 
describes a Bogoliubov sound mode. The second one, 
if present, describes amplitude fluctuations and is associ- 
ated with a Higgs mode. In superfluids, sound excitations 
are gapless while the Higgs mode, if present, is gaped but 
the gap may go to zero under special circumstances such 
as an emergent particle-hole symmetry and Lorentz in- 
variance. This is what happens in the vicinity of the 
superfiuid (SF) to Mott insulator (MI) quantum critical 
point (QCP) of the Bose-Hubbard model when the phase 
transition is crossed at constant density. 

Mean-field theory predicts a stable Higgs particle. In 
(3+1) dimensions, where the QCP is a Gaussian fixed 
point (with logarithmic UV corrections), there is com- 
pelling experimental evidence for the existence of a Higgs 
mode, most beautifully illustrated for the TICUCI3 com- 
pound lj. In (2+1) dimensions, where scaling theory is 
expected to apply, the massive Higgs particle is strongly 
coupled to sound modes and it was argued for a long time, 
on the basis of a 1/iV-expansion to leading order (N = 2 
corresponds to our case), that it cannot survive near crit- 
icality (2j-[5|. Moreover, since the longitudinal suscepti- 
bility diagram has an IR divergence going as a;" 1 , it may 
well dominate any possible Higgs peak. However, it was 
recently emphasized that the type of the probe is impor- 
tant [SHE] : for scalar susceptibility (i.e., the correlation 
function of IV'I 2 ) the spectral function S(uj) vanishes as 
uj 3 at low-frequencies [21 [6] , and this offers better condi- 
tions for revealing the Higgs peak. In the scaling limit 
the theory predicts that S(lu) in the SF phase takes the 
form 

S SF (w) oc A 3 - 2 /^ sf (^) , (1) 



where A is the MI gap for the same amount of detuning 
from the QCP, and v — 0.6717 is the correlation length 
exponent for the U(f) = 0(2) universality in (2 + 1) di- 
mensions [HI HD] ■ The universal function $gp(a;) starts as 
^Sf(^ — > 0) oc x 3 (with small integrated spectral weight 
that may be hard to observe), and saturates to a quasi 
plateau <!>sf(^ > 1) oc x 3 ~ 2 I v « x om25 . The Higgs reso- 
nance (at x ~ 1) can be seen right before the incoherent 
quantum critical continuum with weak w-dependence. 

We are not aware of solid state studies of the Higgs 
mode in two-dimensional (2D) superfluids near the QCP. 
Recently, the cold atom experiment [TT], where a 2D 
Bose-Hubbard system was gently 'shaken' by modulat- 
ing the lattice laser intensity and probed by in-situ single 
site density measurements, saw a broad spectral response 
whose onset softened on approach to the QCP, in line 
with the scaling law ([I]), and no Higgs resonance. This 
outcome can be explained by tight confinement, finite 
temperature, and detuning from the QCP, as shown by 
Quantum Monte Carlo (MC) simulations P2j performed 
for the experimental setup 'as is' in the spirit of the quan- 
tum simulation paradigm [12) . On the other hand, simu- 
lations for the homogeneous Bose-Hubbard model (below 
J, U , and fi stand for the tunneling amplitude, on-site 
interaction, and chemical potential, respectively; in what 
follows energy and frequency are measured in units of J) 

H = - 3 E b h + \ E - 1) - m E ( 2 ) 

<ij> i i 

in the vicinity of the SF-MI point featuring emergent 
particle-hole symmetry and Lorentz-invariance |14j un- 
ambiguously revealed a well-defined Higgs resonance 
which becomes more pronounced on approach to the 
QCP [13]. However, its universal properties, i.e., the pre- 
cise structure of $sf(20j were not answered in Ref. |13| . 
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The Higgs mode is not discussed in the MI phase be- 
cause at the mean-field level there is no Mexican hat po- 
tential to begin with! Likewise, no resonance is expected 
in the normal quantum critical liquid (NL), i.e., at finite 
temperature for critical parameters (U, /x) = (U c ,[jL c ). In 
contradiction with expectations, simulations reveal a res- 
onance in the MI phase right after the gap threshold [13] • 
It was speculated that finite-energy probes are primarily 
sensitive to local correlations at length scales where MI 
and SF are indistinguishable. The universality of the MI 
response was likewise never clarified. 

In their most recent calculation, Sachdev and Podol- 
sky |15j found that including next-order corrections in a 
1 /N expansion in the scaling limit radically changes pre- 
vious conclusions in that S(u>) does contain an oscillatory 
component, in line with MC simulations. However, the 
precise shape of the $gp(x) function could not be estab- 
lished due to artifacts of the approximation. 

In this Letter, we aim to determine the universal scal- 
ing spectral functions when approaching the QCP from 
the SF, MI, and NL phases. We rely on the worm algo- 
rithm [T5HTB] in the path integral representation to per- 
form the required large-scale simulations. By collapsing 
spectral functions evaluated along the trajectories spec- 
ified by the dashed lines in Fig [T] we extract universal 
features for all three phases. They are summarized in 
Fig. [2] which is our main result. Surprisingly, all of them 
include a universal resonance peak (relatively sharp in 
SF and MI phases), followed by a broad secondary peak 
(in SF and MI phases only) before merging with the in- 
coherent critical quasi-plateau (the plateau value is the 
same in all cases, as expected). Our results are in agree- 
ment with scaling theory, and firmly establish that the 
damped Higgs mode survives in all three phases. 

The phase diagram of the 2D Bose-Hubbard model, 
shown in Fig. [I] is known with high accuracy |19H21j 
both at zero and finite temperature. The QCP is located 
at U c = 16.7424(1), fi c = 6.21(2). When the system is 
slightly detuned from the QCP, either by changing the 
chemical potential or the interaction strength, we define 
the corresponding characteristic energy scale A using the 
energy gap in the MI phase, E gap (g), by the rule illus- 
trated in Fig. [lj For positive g = (U — U c ) / J it is half the 
gap, A(g > 0) = E sap (g)/2, where E gap = - fj,i~> is 
deduced from the upper and lower critical chemical po- 
tentials for a given g. For g < along the trajectory i in 
the SF phase it is A(g < 0) = E sap (-g)/2. For U = U c 
and negative g^ = (fi — fi c )/J along the trajectory ii 
in the SF phase we first find g such that fii \g) = fi 
and then define A(g fl ) = CE gap (g)/2 where the constant 
C = 1.2 (see below) is fixed by demanding that the uni- 
versal function is the same along both SF trajectories. 
Note that E gap (g) in the thermodynamic limit can be 
determined accurately from the imaginary time Green 
function data 20j and finite-size scaling analysis. 
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FIG. 1. (Color online) Ground state phase diagram of the 
Bose-Hubbard model in the vicinity of the QCP marked by a 
large (blue) dot (based on Ref. [20] data). The (blue) dashed 
curves specify trajectories in parameter space used to detune 
the system away from the QCP (trajectories i and iii corre- 
spond to unity filling factor n = 1, trajectory ii has constant 
interaction strength). The (black) lines with arrows explain 
how the characteristic energy scale A is obtained for these 
parameters (see text). The inset shows the phase diagram at 
finite temperature, and the trajectory taken in the NL phase. 
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FIG. 2. Universal spectral functions for scalar response in the 
superfluid, Mott insulator, and normal liquid phases. For SF, 
the Higgs peak is at iv H /A = 3.3(8); for MI, tv H /A = 3.2(8); 
and for NL, ojh/T = 6(1). There is a secondary peak around 
uj/ A « 15 in the SF and MI phases, and all responses reach a 
quasi-plateau at the same height 0.6(1) at higher frequencies. 
The error bars on $sf,mi come from the spread of collapsed 
curves, while the ones on "^nl are based on the variance of 
the analytical continuation results |13| . 

To study the scalar response, we add a small uniform 
modulation term to the Hamiltonian Eq.([2j), 

A T 

5H{t) = -5Jcos{ut) bfy = — K(t) , (3) 

<ij> 

where SJ/J <C L The imaginary time correlation func- 
tion for kinetic energy, x( T ) = (K{ T )K(fy) ~ (^) 2 , is 
related to S(ui) through the spectral integral with the 
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finite-temperature kernel, N(t, oj) 

f+OO 



X(r) 



N{t,uj)S(uj) 



-w(1/T-t). 



(4) 



We employ the same protocol of collecting and analyzing 
data as in Ref. [T3]. More specifically, in the MC simu- 
lation we collect statistics for the correlation function at 
Matsubara frequencies uj n — 2ttTti with integer n 



X (Wn) = (K(T)K(0)) t 



(K) 



(5) 



which is related to x( r ) by a Fourier transform. In the 
path integral representation, x(iw n ) has a direct unbi- 
ased estimator, | e luJ " Tk | 2 , where the sum runs over 
all hopping transitions in a given configuration. Once 
X(t) is recovered from x(iw„), the analytical continuation 
methods described in Ref. [T3] are applied to extract the 
spectral function S(uj). A discussion on the reproducibil- 
ity of the analytically continued results for this type of 
problems can also be found in Ref. [13]. 

We consider system sizes significantly larger than the 
correlation length by a factor of at least four to ensure 
that our results are effectively in the thermodynamic 
limit. Furthermore, for the SF and MI phases, we set 
the temperature T = 1/(3 to be much smaller than the 
characteristic Higgs energy, so that no details in the rel- 
evant energy part of spectral function are missed. 

We consider two paths in the SF phase to approach the 
QCP: by increasing the interaction U — > U c at unity fill- 
ing factor n — 1 (trajectory i perpendicular to the phase 
boundary in Fig[T]), and by increasing /i — > /i c while keep- 
ing U — U c constant (trajectory ii tangential to the phase 
boundary in Fig[T]). We start with trajectory i by consid- 
ering three parameter sets for (\g\, L, /?): (0.2424, 20, 10), 
(0.0924,40,20, and (0.0462,80,40). The prime data are 
shown in the inset of Fig. [3] After rescaling results ac- 
cording to Eq. ([I]) we observe data collapse shown in the 
main panel of Fig(3] This defines the universal spectral 
function in the supcrfluid phase, &sf- 

When approaching the QCP along trajectory ii, 
with (\ gii \,L,P) = (0.40,25,15), (0.30,30,15), and 
(0.20, 40, 20) we observe a similar data collapse and arrive 
at the same universal function $sf> see Fig. [4] The fi- 
nal match is possible only when the characteristic energy 
scale A(g^) = CA{g{g^)) involves a factor of C = 1.2. 

The universal spectral function $sf has three distinct 
features: a) A pronounced peak at luh/A « 3.3, which 
is associated with the Higgs resonance. Since the peak's 
width 7/ A m 1 is comparable to its energy, the Higgs 
mode is strongly damped. It can be identified as a 
well-defined particle only in a moving reference frame; 
b) A minimum and another broad maximum between 
w/A 6 [5,25] which may originate from multi-Higgs ex- 
citations P3] ; c) The onset of the quantum critical quasi- 
plateau, in agreement with the scaling hypothesis ([I]), 
starting at lj/A « 25. These features are captured by an 
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FIG. 3. (Color online) Collapse of spectral functions for dif- 
ferent values of U along trajectory i in the SF phase, labeled 
by the detuning g = (U — U c )/J- Inset: original data for 
Ssf(^). 
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FIG. 4. (Color online) Collapse of spectral functions for dif- 
ferent jj, along trajectory ii in the SF phase, labeled by the 
detuning — (/i — g c )/J- Inset: original data for Ssf(u;). 
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We only claim that a plateau is consistent with our imag- 
inary time data and emerges from the analytic continua- 
tion procedure which penalizes gradients in the spectral 
function; i.e., other analytic continuation methods may 
produce a different (oscillating) behavior in the same fre- 
quency range that is also consistent with the imaginary 
time data. 

We now switch to the MI phase, where we approach the 
QCP along trajectory iii in Fig[T] The scaling hypothesis 
for the spectral function has a similar structure to the one 
in Eq. 0, 

Sm(u) cx A 3 " 2 /^ 



'Ml(x) • 



'A' 



(7) 



The low-energy behavior of $mi starts with the thresh- 
old singularity at the particle-hole gap value, &mi(x) ~ 
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l/log 2 (4/(ir-2))6>(2;-2), see Ref. [15]. At high frequen- 
cies $mi(^ 3> 1) has to approach the universal quantum 
critical quasi-plateau (same as in the SF phase). Our 
results for the spectral functions at g = 0.2576 (with 
L = 20,/3 = 10) and g = 0.1276 (with L = 40, /3 = 20) 
are presented in FigjHJ The universal scaling spectral 
function shows an energy gap (this is also fully pro- 
nounced in the imaginary time data). The left side of 
the first peak is much steeper than in the SF phase, in 
agreement with the theoretical prediction for the thresh- 
old singularity. 
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FIG. 5. (Color online) Collapse of the spectral functions for 
different U along trajectory iii in the MI phase, labeled by the 
detuning g = (U — U c )/J. Inset: original data for Smi(w)- 

The universal spectral function in the MI is remark- 
ably similar to its SF counterpart featuring a sharp res- 
onance peak, and thus we interpreted it as a Higgs reso- 
nance. This observation is rather counter-intuitive given 
the standard picture of the Higgs mode in terms of os- 
cillations of the order parameter climbing the hill in the 
Mexican hat potential which is manifestly absent in the 
MI at distances exceeding the correlation length £. Our 
understanding is based on the picture that finite energy 
excitations probe system correlations predominantly in a 
finite space-time volume. The Higgs resonance in the MI 
phase is an indication that for small g and large corre- 
lation length £ oc g~ v the effective action F(^)ji for the 
order parameter coarse-grained over length-scales R <C £ 
still features a Mexican hat potential with strongly renor- 
malized properties. This potential defined under the cor- 
relation length ultimately gives rise to the damped Higgs 
resonance at frequencies exceeding the MI gap. We note 
that since the MI and SF are separated by a critical line 
their scaling functions $mi and 'I'sf remain fundamen- 
tally different at energies smaller than u)jj- 

If the picture of the renormalizcd Mexican hat poten- 
tial is true on the MI side of the transition, it may ex- 
plain why the Higgs resonance was observed at temper- 
atures up to 2 times the Kosterlitz-Thouless transition 
on the superfluid side |13j . even for relatively large de- 



tuning. Furthermore, similar considerations would apply 
to the NL phase at sufficiently low, but finite temper- 
ature T < J (at g = g^ = 0, the superfluid transition 
temperature is zero). In this quantum critical region, 
temperature determines the characteristic energy scale, 
thus 5nl(w) oc T 3 ~ 2 /"$nl(u/T), and all excitations are 
strongly damped. Simulations performed at T/J = 0.5 
on the trajectory iv in the inset of Fig. [I] indeed find 
a peak at low energies before the critical quasi-plateau, 
see Fig.[2j but it is much less pronounced and the oscilla- 
tory component (second peak) is lost. Unfortunately, nu- 
merical complexity did not allow us to verify the scaling 
law directly by collapsing simulations at lower tempera- 
tures and bigger system sizes; we claim the universality 
of ( I?nl(2 ; ) because the value of the plateau is the same 
as in the SF and MI phases. 

In conclusion, we have constructed the universal spec- 
tral functions $ for the kinetic energy correlation func- 
tion for all three phases in the vicinity of the interaction 
driven QCP of the 2D Bose-Hubbard Model. Although 
the nature of excitations in these phases is fundamentally 
different at low temperature, their ^-functions all feature 
a resonance peak which in the SF and MI phases is fol- 
lowed by a broad second peak and evolve then to a quasi- 
platform at higher energy in agreement with scaling pre- 
dictions. In the SF phase, the first peak is interpreted as 
a universal Higgs resonance, which comes from a damped 
Higgs mode. In the MI and NL phase, the existence of 
a sharp peak is unexpected and we interpreted it as a 
Higgs resonance originating from amplitude oscillations 
at mesoscopic length scales. Their origin requires further 
careful theoretical investigation. Experimental verifica- 
tion with cold gases requires flatter traps and lower tem- 
peratures and is accessible within current technology. It 
would signify a new hallmark, going beyond the previous 
studies of criticality near Gaussian fixed points. 
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